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I.  INTRODUCTION 


In  applying  the  Uniform  Geometrical  Theory  of  Diffraction  (UTD)  to 
antenna  radiation  problem  involving  curved  surfaces,  a  major  task  is  to 
determine  the  final  diffraction  point  and  the  geodesic  path  on  the 
curved  surface.  For  the  antennas  mounted  on  the  fuselage  of  an 
aircraft,  the  fuselage  can  be  modeled  as  an  ellipsoid  in  the  UTD 
analysis.  Geodesic  paths  on  an  ellipsoid  have  been  studied  in  detail  in 
References  [2,3]  using  an  elliptic  cone  perturbation  method  which  is 
very  efficient. 

Using  this  perturbation  method  and  another  numerical  technique, 
which  will  be  given  in  this  report,  the  radiation  patterns  for  ellipsoid 
mounted  antennas  is  efficiently  obtained.  The  theoretical  UTD  concept 
to  calculate  actual  radiation  fields  is  given  in  References  [1,21. 

II.  NUMERICAL  TECHNIQUE 
A.  INTRODUCTION 

The  ellipsoid  simulated  by  a  perturbed  elliptic  cone  model  is 
examined  here.  Since  the  elliptic  cone  is  a  developable  surface, 
geodesics  can  be  easily  obtained  [2,3].  Given  a  radiation  direction 
(f)t»  >t).  one  can  find  the  final  diffraction  point  (nQ,  pq)  by  following 
the  geodesic  path,  step  by  step,  until  the  geodesic  tangent  coincides 


1 


with  the  radiation  direction  4>t).  This  is  a  rather  tedius  and  time 

consuming  process  if  applied  for  each  new  radiation  direction. 
Considering  a  new  radiation  direction,  which  does  not  deviate  greatly 
from  the  previous  direction,  one  should  be  able  to  develop  a  solution 
which  uses  the  properties  of  the  surface  and  the  previous  geodesic  path 
to  find  the  new  diffraction  point.  Such  an  approach  is  attempted  here 
to  make  this  solution  as  efficient  as  possible. 

Since  the  field  decays  exponentially  along  the  ray  path  on  the 
surface,  it  is  assumed  that  only  one  or  possibly  two  dominant  rays  exist 
in  the  problems  treated.  One  is  referred  to  References  [2,31  for  more 
details  on  this  topic. 

R.  NUMERICAL  APPROACH  FOR  PATTERN  CALCULATION 

Assuming  the  diffraction  point  is  located  at  0  (a  cos  ve  cos  vr, 
b  cos  ve  sin  vr,  c  sin  ve)  and  the  field  point  at  P  (R^  sin  \  cos  <i>t» 

Rt  sin  ">£  sin  Rt  cos  "'t),  then  at  the  diffraction  point  0  the 
radiation  direction  (Oj.,  should  coincide  with  the  geodesic  tangent  t 
as  shown  in  Figure  1.  Thus, 

t  =  xt/yyzt2 
=  t,  cos  R  +  t  sin  ft 


9 


sin  9.  cos  -  iL_  cos  v  cos  v„ 

*  *  ",  e  r 


sin  o  sin  <>.  -  2_  cos  v  cos  v 


t  vt 


e  r 


cos  9.  -  <L_  sin  v 

t  ,  • 


Note  that 


.2  =  f 


sine,  cos <}>.-,<*  cosv  cosv  ^  +  ^sinB  sin<j>. 
t  t  et  e  r  t  t 


2  2 
-  cosv  sinv  r  +  ' cos 9  -  £_  sinv  1 ' 

e  r  t  ^  e 


=  1  -  2rsin9  cosv  (3  cos<!>.  cosv  +  sin<>.sinv  ) 
tetfT  tr^r  tr 


+  S—  cosS.sinv  1  +  Tcos^v  'll  cos^v  +  sin^v  1  +  si  sin^v  1 
RTte  ep-  r  r?  r  r?  e 

L  t  t  t 


,  and 


ti  =  te  x  n 


*  2222  *  2222 
-xasinvpfb  sin  ve+c 'cos  ve i  +  ybcosvp^a  si n  "ve+c  cos  ve] 

ra2b2sin2v  +c2cos2y  (a2sin2v  +b2cos2v 

e  e v  r  r 


A  2  2 

+zc(b  -a  lsinvr  cosvr  sinve  cosve 

:c2cos2y  +sin2y  fa2cos2v  +b2sin2v  ]U/^ 
e  e  p  p 


t 
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’1 


2  2  2 

bcosvr  cosB(a  sin  ve  +  c  cosve) 

2vr+b2cos^vp) • 

•  [ c2cos2ve+si  n2vg (a2cos2vr+b2si  n2vr) 1 l/2 
bsinve  sinvr  si  n  B 

fc2cos2v  +sin2v  (a2cos2v  +b^sin2v  )U/- 
e  e  r  r 

sin  sin:  -  J?  cosv„  sinv 

t  t  Rt  e  r  .  (2) 

_ ^ 


2  2 

c(b  -a  )sinvr  cosvr  sinve  cosve  cosB 
a2b2si n2vp  +-c^cos2ve  (a2si  n2vr+b2cos2vr)  • 

•  ' c2cos2vg+si n  2ve(a2cos2vr+b2sin2vr) l1/2 
c  cosve  si  nB 

' c 2 c o s  2 v  +s i n  2V  (a 2  cos  ^ v  +b2s in 2  v  ) 1 1 / 2 
e  e  r  r 


cosB  -  n  sinv 
t  dt  e 

- n - 


(3) 


When  the  source  is  located  at  the  off-mid  section  (z  *  0  in  Figure 
1),  the  ellipsoid  is  modeled  by  a  perturbed  elliptic  cone.  The 
associated  unfolded  surface  is  shown  in  Figure  2(b).  If  y  and  B  denote 
the  angle  between  t  and  t.[  at  Q'  and  0,  respectively,  it  is  seen  that 
x.  With  some  manipulation,  one  can  show  that  the  perturbed 
geodesic  path  can  be  expressed  as  follows: 
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re  cos  (y -n)  =  rs  cos  y 


(4) 


where 


«(Vp) 


r  fa^b?  +  Z?  cot4 V  (a^sin^v'  +  b^cos^v'lU/? 

•  s  s  _s  es  s  r  s  r  dV 

Vps  a^cos^v'  +  b^sin^v'  +  Z^cot4\/  r 

s  r  s  r  s  e* 


a„  =  acosV  ,  b  =  bcosV  ,  Z„  =  csirW  , 
s  e„  s  e„  s  e„  ’ 

s  s  s 


/  ^  2.,  .2  .  2tJ  .  -.2  .4..  1/2 

rs  =  (ascos  Vr  +  bssin  Vr  +  Zscot  Ve  ) 


7  7  7  7  7  A  1  /? 

re  =  (ascos^Vr  +  b$sin  Vr  +  Zgcot  v6s )  -  $e 


and 


C  ,Ve  r  2  2  ‘  ,  2  2.,  k2  .  2W  .  .  2  'i1/*' 

=  f  c  cos  V  +  (a  cos  V  +  b  sin  V  )sin  V  !  dV 
e  Jw  1  e  r  r  e  e 


Now,  [tx  (-b  sin  Vr)  +  ty  a  cos  Vr]  yields 


7  7  7  7  7  7  7  1/2 

ab  cos(y-o)  [c  cos  Ve  +  sincVe  (a  cos  Vr  +  b  sin  Vr)l 

[a^h^sin^V^  +  c^cos^Ve  (a^sin^V^  +  b^cos^V^)]!^ 


sin  (acosVr  s i n <>t  -  bsinVr  cos^) 
- D - 


Next,  |[txbcosVr  +  tyasi nVrlccosVe  +  tzabsinVe'  yields 


(5) 


a 


absinVe  cos9t  +  csin9t  cosVe  (asin^sinVp  +  bcos^cosVp) 

- D - 


_  abc  =  o  (6) 

w 

Three  functions  can,  then,  be  constructed  as  follows  from  Equations  (4)- 

(6): 

F (Ve,Vr, y)  =  recos(y-n)  -  rscosy  -  0  (7) 

^t *  ft »  Ve,  Vr,  y) 

=  Dabcos  ( y-oi)[c2cos^Ve  +  sin2ve(a^cos^Vr  +  b^sin^Vp)!  1/2 

-  sin9t  (asin^  cosVr  -  bcos^  sinVr)  • 

•  [a2b2sin2Ve  +  c^cos^Vg  (a^sin^Vp  +  b^cos^Vr)Tl/2 

=  0  .  (8) 

Further,  one  finds  that 

8(Rt*  ®t*  ‘ft *  ^e»  ^r) 

=  absinVe  cosOt  +  csin9t  cosVe  (asin^  sinVr  +  bcos«t  cosVr) 

-  abc  =  o  .  (g) 


q 


Provided  that  one  has  obtained  a  diffraction  point  (Ve>  Vr)  for  a 


receiver  location  (R^,  'It)*  a  numerical  technique  can  now  be 

developed  from  Equations  (7),  (8),  and  (9)  to  solve  for  (Ve  +  A Ve,  Vr  + 
AVr)  associated  with  a  new  receiver  location  (R^  +  aR^,  +  ao^,  ^  + 
a>(.).  Assuming  that  the  ith  set  of  (R^,  />t»  ve»  Vr)  is  first  known 

to  satisfy  F-j  =  H-j  =  G-j  =  0,  or  at  least  approximately  so,  the  next  set 
(Rt  +  \Rt,  >  +  A‘\,  ^  +  A t>t ,  Ve  +  a Ve,  Vr  +  AVr)  is  obtained  by 
enforcing  F j  =  H-j+j  =  Gj  =  0,  such  that 

F .  .  =  F .  +  Fw  AV  +  Fw  AV  +  F  Ay  =  0  , 

i+l  i  V  e  V  r  y 
e  r 


G =  G.  +  Gu  AV  +  G..  aV  +  G  Ay 
i+l  l  V  e  Mr  y 
e  r 


+  G  (  \  ’  +  G  s  A->.  +  Gp  ARt  =  0 

t  c  i  Kt 


and 


H.  .  =  H.  +  H„  AV  +  AV  +  H  A  A 
i+l  i  Me  Mr  \  t 
e  r  t 


+  H  ,  A  A. 
't  1 


aR. 


=  0 


In  matrix  form,  it  is  given  by 


FV  FV  Fy 

AV 

-Fi 

e  r 

e 

Gw  Gw  Gv 

A  V 

-G,  -  Gn  An  -  Gi,  A  a  -  Gr  Ar 

e  r 

r 

t  t  t  pt  t  t 

Hw  Hw  0 

Ay 

-Hi  -  Hn  An  -  A*  -  Hr  Ar 

e  r 

t  t  t  t  t  t 

(10) 
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2  2  2 
Dab  cos(y-a)(b  -a  )sinVr  cosVr  sin'Ve 

r  [c^cos^V  +(a^cos^V  +b2sin^V  )sin2v  ]l/2 
e  r  r  e 

+  Dab[c2cos2Ve+(a2cos 2Vr+b2si n2Vr )si n2Vel^ ^2  sin(Y-'0  . 

[a^b^+Z^cot4Ves(a^sin2Vr+b^cos2Vr)l1/2 

a^cos^V  +b2sin2v  +z2cot4V 
s  r  s  r  s  es 

2  2  2  2 

si  n  (asi  n  ^cosVr-bcoS‘t>tsinVr)  (b ‘-a  )si  nVrcosVrc  cos  Ve 

[(a2sin2V  +b2cos2V  )c2cos  +a^b^sin^V 
r  r  e  e 

+  sin-\(asin  ^sinVr+bcosi>tcosVr)r(a^sin^Vr+b2cos^Vr)  • 

•  c^cos^Vg+a^b^sin^Vg]1/2 

+  3bcos  ( y-n)  c2cos2V  +(a2cos2V  +b2sin2V  )sin2V  ]^2  • 

n  e  r  r  e 

.  sinV  cos V  cos2V  -sin9  cosV  (&  si n <{>.  cosV 

r  r  e  t  vt  r 

t  t 

-  a  cos  s i  n V  }  1 
*t  1  r 


(i ,  =  -Dabsio(Y- T)[c^cos^Vg+(a^cos^Vr+b2sin2vr)sin2vell/2 


G9t  =  -  cose^asin^cosI^-bcos^sinV^)  • 


•  [(a2sin2Vr+b2cos2Vr)c2cos2Ve+a2b2sin2Ve]1/2 


+  abcos(Y-a)  [c2cos2Ve+(a2cos2Vr+b2sin2Vr)sin2Vel1/2 


{  sin0tsinVe-cos6tcosVe(3_  cos<)>tcosVr+  b__  sin<t>tsinVr)i 


=  -  sin9t(acos4>t.cosVr.+bsin<ttsinVr)  • 

•  [(a2sin2Vr+b2cos2Vr)c2cos2Ve+a2b2sin2Ve]1,/2 

+  abcos(v-a)  C c2r.os2V  +f  a2cos2V  +b2sin2V  )sin2V  ]1//2  • 

IT  e  r  r  e 


si  no.  cosV 
t  e 


(£ 


'*7 


sin<>.cos V  - 
t  r 


b 

*7 


cos^sinV^)' 


cos(y-a)[c2cos2Ve+(a2cos2Vr+b2sin2Vr)sin2Ve]^2  • 


.  {[sine  cosV  (a  cos$tcosV  +b  sin4>tsinV  )  +  c  cose  sinV  ]/„ 

vc  ur  vr 


-  [cos2Ve(a2cos2Vr+b2sin2Vr)+c2sin2Ve]/R3' 
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Hy  -  ah  cosV  cos 9  -  c  sinV  sine  (asin^.sinV  +  bcos^.cosV  ) 

6  sl  ci  t  r  i  r 

Hyr  -  c  cosVg  sinS^.  (asi n^cosVr  -  bcos^sinV^ 

Hy  =  0 

H  s  -ab  sinV  sino.  +  c  cosV  cos 9.  (asi n <j>.  si nV  +  bcos^.cosV  ) 

j  G  L  Cl  l  I  l  r 

H*,  =  c  cosVesinet  (a  cos<{>tsinVr  -  b  si  n^cosV^ 

and 

Hr  -  abc 
t  ^r 
t 

It  is  seen  that  one  can  solve  for  (a Ve,  \Vr,  Ay),  for  a  known  ( aR^ . 

\ ^ ) ,  using  Equation  (10).  To  obtain  a  diffraction  point  (Ve,  Vr)  for  a 
given  receiver  location  (R^ ,  et,  <^),  one  can  always  assume  the  first 
diffraction  point  is  at  the  source  (Ve,  Vr)  =  (ves,  Vr$)  with  the 
radiation  direction  (9^,  =  ^.)  for  the  positive  ray  (in  Y  direction) 

or  ("y,  Rj.  =  3_n)  for  negative  ray  (in  -Y  direction),  and  gradually  add 
the  increments  (aR^,  A9^,  a^)  until  the  final  radiation  direction 
(">t»  ?t)  is  reached  as  shown  in  Figure  3.  More  detail  on  this  topic  is 
provided  in  Reference  [4],  One  need  not  start  out  from  the  source 
everytime,  but  obtains  the  new  diffraction  point  directly  from  Equation 
(10),  provided  that  the  new  receiver  location  does  not  deviate  greatly 
from  the  previous  direction. 

After  the  geodesic  path  is  determined,  various  other  parameters 
associated  with  actual  field  calculation  must  be  found.  The  Fock 
parameter  r  was  calculated  in  Reference  [2]  as  follows: 


I 


=  rscosy 


Vr 

[  ^  ,)g  1  1  /  3  _  1 _ (j;, 

Vps  °g  ^  cos ^  {  y  -  x )  d\T  r 


where 


du 


O  O  9  4  ^  9  i  9  9  | 

[a  h  +Z  cot  ve  (a  sin  V  ♦  b  cos'V  )"| 
sss  s  s  r  s  r 

a^cos^v'  +  n^sin^V  +  "z^cot 

s  r  s  r  s  es 


1/? 


or 


•vc 


1  koqU/3 


sin ( y-o) 


dS 

dV' 


e  d  V 


Note  that  o 

curvatures . 


_ _  1 _ and  k.  and  k„  are  two  principal 

9  k^ cos ^ 3  +  k^sin^Z 


Next,  the  ray  divergence  factor  ft is  defined  as  the  change 

/  dp(O’) 

in  the  separation  of  adjacent  surface  rays  as  shown  in  Figure  4.  Since 
the  ellipsoid  simulating  the  aircraft  fuselage  will  be  long  and  slender, 
it  is  assumed  that  the  ray  divergence  factor  is  unity  in  the  analysis. 
This  completes  the  elliptic  cone  perturbation  solution  for  the 


antenna  mounted  on  the  off-mid  section  of  an  ellipsoid. 


Figure  4.  Illustration  of  the  divergence  factor  (/d\|*0/d4>) 
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III.  RESULTS 


The  solutions  presented  in  the  previous  chapter  are  employed  to 
compute  the  near  field  radiation  patterns  for  short  monopoles  or  slots 
mounted  on  an  ellipsoid. 

To  examine  different  conical  pattern  cuts,  a  cartesian  coordinate 
system  (x',y',z‘)  originally  defining  the  ellipsoid  geometry  is  now 
rotated  into  a  new  system  (x,y,z)  as  shown  in  Figure  5.  Note  that  the 
new  cartesian  coordinates  are  found  by  first  rotating  about  the  z'-axis 
a  angle  pc  and  then  about  the  y-3xis  a  angle  hc.  The  pattern  is,  then, 
taken  in  the  (x,y,z)  coordinate  system  with  Op  fixed  and  -v  varied. 

To  show  the  validity  of  the  elliptic  cone  perturbation  solution, 
some  typical  sources,  i.e.,  short  monopole,  axial  slot  and 
ci rcumferential  slot,  and  various  source  locations  are  chosen  as  shown 
in  Figure  6. 

For  each  case  the  following  typical  radiation  patterns  are 


obtained 

• 

a ) 

'C 

=  0°, 

*c  ; 

=  90°, 

!P  = 

90° 

(roll  plane 

pattern ) 

b) 

\ 

=  30°, 

% 

=  90°, 

°P 

=  90° 

c) 

=  60°, 

% 

=  90°, 

=  90 c 

d) 

}c 

=  90% 

% 

=  90°, 

’3P 

=  90° 

1  (elevation 

plane  pattern) 

e) 

=  90°, 

':’c 

=  o°. 

.>  - 

P  " 

90° 

(azimuth  pic 

me  pattern). 

The  radiation  patterns  obtained  by  the  ellipsoid  program,  which 


uses  an  ellipsoid  to  simulate  the  aircraft  fuselage,  are  compared  to 
those  obtained  using  the  spheriod  solution  [51  in  each  case. 


near  field 
receiver 

LOCATION 


Definition  of  pattern  axis 


It  is  noted  that  the  geodesic  tracing  method  of  the  ellipsoid 
program  for  the  side  mounted  antennas  (Figures  9,  10,  11,  14,  15,  18, 

19)  is  different  from  that  of  the  spheroid  program  because  the  ellipsoid 
is  not  a  surface  of  revolution. 

The  exact  agreement  between  the  results  of  the  ellipsoid  program 
and  the  spheroid  program  as  shown  in  Figures  7-19  gives  one  confidence 
about  the  validity  of  the  elliptic  cone  technique. 

Next,  the  ellipsoid  program  is  employed  to  calculate  the  radiation 
patterns  due  to  antennas  mounted  on  an  ellipsoid  surface.  The  typical 
ellipsoid  geometry  (2a  x  4a  x  10 \ )  is  chosen  and  examined  for  various 
sources  and  source  locations. 

The  cone  boundary  shown  in  Figure  32  is  used  in  determining  whether 
one  or  two  rays  are  used  in  the  solution.  Note  that  ,i2  is  defined 
automati cal ly  by  determining  the  caustic  angle  in  the  elevation  pattern 
(■’c)  and  adding  a  few  additional  degrees  to  that  value,  i.e.,  *1?  = 

3C  +  A  ’  where  2°  s  \,t  s  10°.  One  would  expect  to  observe  slight 
discontinuities  somewhere,  because  various  numbers  of  rays  are  included 
in  different  regions. 

IV.  CONCLUSIONS 

The  object  of  this  study  has  been  to  develop  an  efficient  numerical 
solution  for  the  high  frequency  radiation  patterns  of  an  ellipsoid- 
mounted  antenna.  The  UTO  is  used  in  this  study  to  calculate  the 
radiation  patterns,  and  the  elliptic  cone  perturbation  method  is  applied 


to  simulate  the  geodesic  paths  on  the  ellipsoid,  which  in  turn  can  be 
used  to  model  an  aircraft  or  missile  fuselage.  For  a  given  radiation 
direction  in  the  shadow  region,  the  geodesic  path  and  final  diffraction 
point  on  the  ellipsoid  can,  then,  be  found  via  an  efficient  numerical 
approach. 

The  exact  agreement  of  the  radiation  patterns  from  two  different 
programs  confirms  that  this  elliptic  cone  perturbation  solution  is  very 
useful  in  predicting  the  high  frequency  radiation  patterns  for  antennas 
mounted  on  the  off-mid  section  of  an  ellipsoid. 

This  numerical  solution  will  be  employed,  along  with  flat  plates  to 
construct  a  general  solution  for  calculating  radiation  patterns  due  to 
airborne  antennas. 


?? 


180° 


(d)0c  =90°  ,  ^>c  =90°  ,  @p  =90° 


Ellipsoid  Program  Spheroid  Program 


Figure  7.  (continued) 


ld)9c  -90°  ,  <^c  =90° ,  dp  =90° 


(e)0c=9O°,^c  =0° ,  0p  =90  ° 


Ellipsoid  Program 


Spheroid  Program 


fiqurp  H.  (continued) 


?f> 


td]0c  =90°  ,  <^>c  =90°  ,  =90° 


(e)  0  c  =  90°,  ^c=O°,0p=9O° 


Ellipsoid  Program 


Spheroid  Program 


Figure  11).  (continued) 


0° 


(b)8c  =30°, 


4>c=9O°,0p=9O° 


i 


/ 


Ellipsoid  Program 


Spheroid  Program 


figure  12.  Comparison  of  radiation  patterns  for  an  axial 
slot  mounted  at  *s  =0°,  Bs  =  60°  on  a  2  \  x 
1 0 x  spheriod. 
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Id)  6C  =90°,<pc  =90°  ,  9p  =90° 


le!  0C=9O°,  <£c=O°,0p  -90° 


Ellipsoid  Program 


Spheroid  Program 


l  i  i|ur.>  1  / .  ( root  i  ntiprl  1 
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Ellipsoid  Program 


Spheroid  Program 


Figure  13.  (cont  i niiprl ) 
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(e)  0c  =  9O°,<£c=O°(0p  =90° 


EUipsoid  Program 


Spheroid  Program 


Figure  15.  (continued) 


Ellipsoid  Program 


Sphere 


Figure  16.  (continued) 


(d)  ec=9o°,  4,^90°,  ep  =90° 


Ellipsoid  Program 


Spheroid  Program 
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Id.  (continued) 
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